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' Context. The remarkable properties of CS 30322-023 became apparent during the course of a high-resolution spectroscopic study of a sample 

of 23 carbon-enhanced, metal-poor (CEMP) stars. 
■ Aims. This sample is studied in order to gain a better understanding of s- and r-process nucleosynthesis at low metallicity, and to investigate 

the role of duplicity. 

. Methods. High-resolution UVES spectra have been obtained, and abundances are derived using 1-D, plane-parallel OSMARCS models under 
' the LTE hypothesis. The derived atmospheric parameters and observed abundances are compared to evolutionary tracks and nucleosynthesis 
^ (— i predictions to infer the evolutionary status of CS 30322-023. 

Q | Results. CS 30322-023 is remarkable in having the lowest surface gravity (logg < -0.3) among the metal-poor stars studied to date. As a 
result of its rather low temperature (4100 K), abundances could be derived for 35 chemical elements; the abundance pattern of CS 30322-023 
is one of the most well-specified of all known extremely metal-poor stars. With [Fe/H]= -3.5, CS 30322-023 is the most metal-poor star to 
exhibit a clear s-process signature, and the most metal-poor "lead star" known. The available evidence indicates that CS 30322-023 is presently 
a thermally-pulsing asymptotic giant branch (TP-AGB) star, with no strong indication of binarity thus far (although a signal of period 192 d 
is clearly present in the radial-velocity data, this is likely due to pulsation of the stellar envelope). Low-mass TP-AGB stars are not expected 
to be exceedingly rare in a magnitude-limited sample such as the HK survey, because their high luminosities make it possible to sample 
them over a very large volume. The strong N overabundance and the low l2 C/ l3 C ratio (4) in this star is typical of the operation of the CN 
cycle. Coupled with a Na overabundance and the absence of a strong C overabundance, this pattern seems to imply that hot-bottom burning 
C$ ' operated in this star, which should then have a mass of at least 2 M Q . However, the luminosity associated with this mass would put the star at 
a distance of about 50 kpc, in the outskirts of the galactic halo, where no recent star formation is expected to have taken place. We explore 
alternative scenarios in which the observed abundance pattern results from some mixing mechanism yet to be identified occurring in a single 
low-metallicity 0.8 M B AGB star, or from pollution by matter from an intermediate-mass AGB companion which has undergone hot-bottom 
burning. We stress, however, that our abundances may be subject to uncertainties due to NLTE or 3D granulation effects which were not taken 
into consideration. 
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1. Introduction 

The star CS 30322-023 first appeared in the lNorris etai]dl999l) 
list of broadband photometry for candidate met al-poo r stars se- 
lected from the HK survey dBeers et al.lll985l Il992l) . Follow- 
up medium-resolution spectroscopy of CS 30322-023 revealed 
that it was likely to be rather C-rich. It was therefore included 
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in our sample of 23 carbon-enhanced metal-poor (CEMP) stars 
observed with the UVES spectrograph on VLT-Kueyen, in or- 
der to look for more lead s tars after the initial discoveries by 
IVanEcketallfcOOll 2003). Lead stars are low-metallicity ob- 
jects where the operation of the s-process nucleosynthesis was 
so efficient that it converted a rela tively large fraction of iron - 
seed nuclei all the way to lead JGorielv & Mowlavil 12000). 
The sample we investigated was collected from three different 
sources: 

- CEMP objects from the HK survey not included in high- 
resolution follow-up studies by other teams; 
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- Stars from the list of faint high-latitu de carbon star s in the 
Hamburg-ESO (HE) survey JChristlieb etalll200lh which 
have strong CN or C2 bands; 

- The catalog ue of low-metallicity stars w ith strong CH 
bands from lsieivvte & Bartkeviciusl dl99d) . 

Prelimina ry an a lyses of this sample have been presented 
bv iMasseron et al] d2005l |2006); a more detailed paper is in 
preparation as part of Masseron's Ph.D. thesis {Observatoire 
de Paris, France). 

The star CS 30322-023 deserves, however, a dedicated 
analysis, as it exhibits several pecularities as compared to other 
metal-poor stars. First, as will be shown in Sect. [6] it is an- 
other example of a lead star. Secondly, despite the presence 
of a definite, albeit moderate, s-process signature, it has al- 
most no C enhancement, the strength of the CN bands be- 
ing due to a strong N overabundance. It is therefore one of 
the few known Nitrogen-Enhanced Metal-Poor (NEMP) stars. 
Thirdly, its most intriguing property is, by far, its very low 
gravity (log g about -0.3), coupled with a rather low temper- 
ature (4100 K; Sect. [3}. Its evolutionary status is thus puzzling 
(Sect. @J; is it a low-mass (~ 0.8 M Q ) star close to the tip of 
the AGB or a rather luminous intermediate-mass star located 
in the outskirts of the galactic halo? Either possibility raises in- 
triguing questions concerning its nature, which are addressed 
in Sects. |5]andQ Among these are: What is the probability of 
catching a low-mass star in an evolutionary stage so close to 
the tip of the AGB? Is the star formation history in the halo 
consistent with late formation of an extremely low-metallicity 
intermediate-mass star in the outskirsts of the Galaxy? Is it a 
star accreted during a merger event with a satellite galaxy? Is it 
a star that evaporated from a globular cluster? Is it a low-mass 
star polluted by an intermediate-mass companion formerly on 
the thermally-pulsing aymptotic giant branch? 

2. Observations and radial-velocity variations 

Several high-resolution spectra of CS 30322-023 have been ob- 
tained at the European Southern Observatory with the UVES 
spectrograph mounted on the VLT-Kueyen telescope, and with 
the FEROS spectrograph mounted on the 2.2m telescope (cour- 
tesy S. Lucatello). The observation log is listed in Tabled The 
first and last spectra have been obtained in the framework of 
other UVES progra mmes (such as the ESO Large Programme 
on "First Stars"; e.g.. lSpite et al.l2o'()6l) . and the corresponding 
radial velocities have been kindly communicated to us by M. 
Spite and J. Johnson. The two spectra exposed for 3600 s on the 
blue arm of UVES (two times 1800 s on the red arm) were used 
to derive the abundances (Sect.[6j. Their signal-to-noise (SNR) 
ratio is around 200/1 at 400 nm. The spectra exposed for 90 s 
were only used to derive radial velocities by cross-correlating 
the observed spectrum with the two high-SNR spectra. The 
adopted heliocentric velocity for a given low-SNR spectrum 
is the average of these two differential values; it is listed in 
Tabled Arcturus and solar templates were used to estimate the 
internal consistency of the radial-velocity measurement, and a 
telluric template spectrum was used to evaluate the zero point 
offset velocity. The uncertainty on the radial velocity is the root 



Table 2. Available photometry for CS 30322-023; m and m 
are the observed and de-reddened magnitudes in the considered 
band, respectively. 



Band 




+ 


m 


Ref. 


U 


14.48 


0.02 


14.19 


(1) 


B 


13.49 


0.03 


13.24 


(1) 


V 


12.21 


0.03 


12.02 


(1) 


./ 


10.050 


0.024 


10.00 


(2) 


H 


9.479 


0.024 


9.45 


(2) 


K 


9.313 


0.021 


9.30 


(2) 


Bol. 






11.7 


(3) 



(1) Norris et al. 1 1999); (2) Cutri et al. 1 2003); (3) This work 

mean square of the template error and the precision of the cen- 
tering of the gaussian fit to the cross-correlation profile. Both 
uncertainty sources, along with their root mean square error, 
are listed in Tabled 

The standard deviation of the observed radial velocities 
amounts to 1.24 km s _1 , which is considerably larger than the 
average uncertainty, amounting to only 0.18 km s _1 . Possible 
causes for these variations are investigated in Sect. 14.2] 

Broadband UBV photometry is provided by Norr is etall 
dl999l) . while near-infrared JHK photometry is available from 
the Two Micron All Sky Survey dCutri et al.ll2003l) . as listed in 
Table |3 Given the large distance to CS 30322-023 (Sect.gJ, 
the redden ing along its line of sigh t across the Galaxy was 
taken from Burstein & Heiles ( 1982) to be Eb-v = 0.06, from 
which it follows that Ay = 0.18 mag . This value is close to the 
0. 16 mag predicted by the model of lArenou et alJ {l992). The 
reddening law A^/Ay from lCohen et alJ (119811) . appropriate for 
carbon stars, has been applied to de redden the observed ma gni- 
tudes. Zero-magnitu de fluxes from Schmi dt-Kaler! fl982) and 
Coh en et al.l ([2003b ) have been used to construct the spectral 
energy distribution. The bolometric flux has then been obtained 
by computing the integral J AFx din A using the trapezoidal 
rule on the UBV and JHK bands. A bolometric correction of 
-0.30 mag is finally obtained, yielding Wboi = 11.7. It should 
be noted that this bolometric correction differs so mewhat f rom 
the v alue of -0.8 mag predicted by Table 6 of Alonso et al. 
dl999l) for giants with [Fe/H]= -3 and r eff = 4100 K. 

3. Atmospheric parameters of CS 30322-023 

The effective temperature of CS 30322-023 was obtained, in 
the usual manner, by requiring that the abundances derived 
from various iron lines exhibit no trend with excitation po- 
tential. Gravity was obtained through ionisation balance, and 
microturbulence was set by requiring the absence of a trend 
between abundances and reduced equivalent widths. A photo- 
metric 7W es timate has been obtained from the calibration of 
lAlonso et all dl998l) . applied on the 2MASS / - H and J - K 
indices, assuming a metallicity of -3.0. 

The derived atmospheric parameters, usin g model atmo- 
spher es belonging to the OSMARCS family ( Gusta fsson et alJ 
2003), are listed in Table[5] The discrepancy between the pho- 
tometric and spectroscopic temperature estimates probably re- 
sults, on one hand, from the fact that the CS 30322-023 param- 
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Table 1. Observation log. The uncertainty on the radial-velocity Vr arising from the template mismatch (e temp ), from the gaussian 
fit of the correlation profile (e/, f ), and their root-mean-square (e tot ), are also listed. 



Date 


HJD 


Exp. time 


Vr 


€tem p 




Etot 


Telescope 


Programme 


(YYYY/MM/DD) 


(-2 450 000) 


(s) 


(km/s) 


(km/s) 


(km/s) 


(km/s) 






2001/11/04 


2218.49762 


3100 


116.53° 




0.08 


0.08 


VLT/UVES 


165.N-0276(A) 


2002/04/24 


2388.83689 


3600 


116.40 


0.04 


0.04 


0.06 


VLT/UVES 


069.D-0063(A) 


2002/05/11 


2405.86677 


3600 


114.59 


0.05 


0.04 


0.06 


VLT/UVES 


069.D-0063(A) 


2004/04/29 


3124.92121 


90 


115.92 


0.04 


0.04 


0.06 


VLT/UVES 


073.D-0193(B) 


2004/06/14 


3170.81670 


90 


115.75 


0.04 


0.03 


0.05 


VLT/UVES 


073.D-0193(B) 


2004/08/17 


3234.64390 


90 


112.81 


0.05 


0.02 


0.05 


VLT/UVES 


074.D-0228(A) 


2004/11/11 


3320.53675 


90 


115.71 


0.05 


0.03 


0.06 


VLT/UVES 


074.D-0228(A) 


2005/06/29 


3550.89234 


900 


115.09* 


0.37 


0.05 


0.37 


2.2m/FEROS 


075 .D-0 120(A) 


2005/07/26 


3577.86729 


900 


114.61* 


0.39 


0.08 


0.40 


2.2m/FEROS 


075.D-0120(A) 


2005/07/27 


3578.81266 


900 


113.69* 


0.29 


0.09 


0.30 


2.2m/FEROS 


075.D-0120(A) 


2005/07/28 


3579.89378 


900 


113.82* 


0.61 


0.08 


0.61 


2.2m/FEROS 


075.D-0120(A) 


2005/08/17 


3599.67995 


2000 


113.35 c 


0.03 


0.05 


0.06 


VLT/UVES 


076.D-0451(A) 



" Data kindly communicated by M. Spite 
* Data kindly communicated by S. Lucatello 
c Data kindly communicated by J. Johnson 



eters fall at the edge of the validity domain of Alonso's calibra- 
tion and, on the other hand, from model errors impacting the 
spectroscopic determination, as we now discuss. Fig.n reveais 
that a temperature T e s = 4300 K, as suggested by the pho- 
tometric calibration, combined with a higher gravity, does not 
yield a satisfactory solution. Fig.|2j which presents the Fe I and 
Fe II abundances derived from various lines for plane-parallel, 
LTE model atmospheres, shows that the best solution is ob- 
tained for T e g = 4100 K and logg = -0.3. The uncertainty 
on the gravity, in the framework of Fe I/Fe II ionisation equi- 
librium under LTE (see, however, the dicussion about NLTE 
effects below), may be estimated by comparing the three pan- 
els of Fig. |3 corresponding to logg = 0.5,0.0 and -0.5. The 
gravity of CS 30322-023 clearly lies between -0.5 and 0.0. The 
value log g = -0.3 + 0.3 has therefore been adopted. 

We have also attempted to employ spherical atmosphere 
models, which are better suited to low-gravity stars than the 
conventionally used plane-parallel models. The Fe I/Fe II ioni- 
sation balance would require spherical models with an effective 
gravity logg < -0.3. We could not, however, converge mod- 
els with such low surface gravities, due to the combination of 
the large atmospheric extension and the scattering character of 
the radiation field at such low metallicities, which is decoupled 
from the local thermal pool of the gas. 

It should be noted that the gravities of metal-poor gi- 
ants derived from LTE iron-line analyses are most proba- 
bly u nderestimated, because of the neglect of non-LTE ef- 
fects dThevenin & Idiard 1 19991 llsraelian et all l200ll |2004; 
iKorn et al.l2003h . We plan to investigate these effects in a forth- 
coming paper with models specifically tailored to CS 30322- 
023. Nevertheless, the conclusion that CS 30322-023 is the 
most evolved star found thus far among CEMP stars remains 
valid, as the comparison presented in Fig.[3]involves stars ana- 
lyzed under LTE as well. In absolute terms, with NLTE correc- 
tions on the gravity of 0.5 dex (as derived by Thevenin & Idiart 
Il999t IKorn e7alll2003l for giant stars with [Fe/H]= -3) or 
even up to 1 .0 dex (from the comparison of LTE and NLTE 
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Fig. 1. Iron abundances derived from all lines measured in 
CS 30322-023 (crosses: Fe I lines; open circles: Fe II lines), 
as a function of the lower excitation potential for the adopted 
model (Lower panel: T e s = 4100 K, logg = -0.3) and for a 
model of higher temperature and gravity (Upper panel: T e ff = 
4300 K, log g = 0.5). The latter model - the best one at 4300 K 
- is seen to be less satisfactory than the adopted 4100 K model. 
In both cases, plane-parallel LTE model atmospheres (with 
[Fe/H] = -3.5 and a microturbulent velocity of 2.2 km s -1 ) 
have been used. 
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Fig. 2. Same as Fig^for different gravities at T e ff = 4100 K. 
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Table 3. Atmospheric parameters of CS 30322-023. In the third column, 'PP' stands for plane-parallel and 'S' for spherical (see 
text for a discussion of the uncertainties). 



7W (K) 

microturbulent vel. (km/s) 



4100 excitation balance (PP-model) 

4100 excitation balance (S-model) 

4350 2MASS photometry 

-0.3 ± 0.3 Fe I/Fe II ionisation balance in PP-model (LTE) 

< -0.3 Fe I/Fe II ionisation balance in S-model (LTE) 

2.2 PP-model 

1.7 S-model 



gravities for CS 2 9498-043 derived by lAoki et alJ l2002dt 
llsraelian et ai1l2004 . CS 30322-023 would remain on the TP- 
AGB (Figs.fJandE). 

4. Evolutionary status of CS 30322-023 

4.1. Position in the Hertzsprung-Russell diagram 

The gravity and effective temperatures derived in Sect, [3] may 
be combined with the apparent bolometric magnitude to de- 
rive the distance of CS 30322-023 as a function of the adopted 
mass, using the relationship: 

logM/M = log g/g Q + log L/L@- 4 log T eff /T eft -, Q (1) 
= -2.25 - 0.4M bol (2) 
= 2 log d(pc)- 8.93, (3) 

where we have adopted T e $ = 4100 K, logg = -0.3 and 
m bo i = 11.7 for CS 30322-023 [denoted set (1) in Tableg). The 
corresponding distances and luminosities are listed in Table |4] 
To illustrate the effects of the uncertainties, values for log L/L e , 
Mboi and d [listed under set (2)] have been derived assuming 
that the gravity could be 1 dex higher (due to NLTE effects) 
and the temperature 250 K higher. 

CS 30322-023 appears to be the star with the low- 
est gravity detected so far among CEMP stars, as shown 
in Fig. w hich collects all such stars from the litera- 
ture toearborn et al.lll986l iRvan et alJll99ll iMc William et al 



1995; [Norris et alJll997albt ISneden et alJll99a iGratton et al 
200 0t iHill et alJl2Q0(A ISneden et alJl2000t IWestin et alJl2 000: 
Norris etalJ l200lt IPreston & Snedenl l200lt IVanEcketaL 
200 ll iTravaglio et alJkOOlt lAoki et alJkOOlt IChristlieb et al. 
2002 iHilLetalJ |2002l iJohnson & Bolte 2002; Aoki et al 



2002 alblcldl 2Q03t ICohen et alJl2003at iLucatello etaiT 2003; 



Sne den et al l2003t IVan Eck et alJ 120031 lAoki et alf 2004; 



Cohen et al. 2004; Honda et al. 2004; Johnson & Bolte 2004; 
Simmerer et all2004 iBarbuv et al 120051: IPlez & CoheiJ2005l 
Frebel et all2005tl200d:lMasseron et all200d) . 



With T eB = 4240 K and logg = 0.3, the star HD 110184 
(which, despite the reported detection of Th, is not con- 
sidered an r-process enhanced star, according to the crite- 
ria of Beer s & Christlie b 2005 ; its [Eu/Fe] ratio is too low) 
studied by iHonda et"al ] 12004 is the second most evolved 
star aft er CS 30322-023. Then comes CS 29498-043 stud- 
ied bv lAoki etalJ J2002di l2004 . IWeiss et alJ (I2004 . and 
llsraelian et all d2004 . with r eff = 4400 K and logg = 0.6 
(Aoki et al. 2002d, a subsequent analysis by Aoki et al., 2004, 
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Fig. 3. Hertzsprung-Russell diagram for CEMP stars collected 
from the literature (see references in the text), revealing that 
CS 30322-023 (asterisk) is the most extreme case known. The 
solid line corresponds to the evolutionary track of a 0.8 M 
model with [Fe/H] = -3.8 (see text). The horizontal bar on 
the evolutionary track marks the tip of the RGB. The arrow 
corresponds to a conservative 250 K uncertainty on r e ff and a 
1 .0 dex uncertainty on the gravity of CS 30322-023 due to the 
use of LTE instead of NLTE model atmospheres (see text). 



yielded instead T eS = 4600 K and l og r = 1.2). Using a NLTE 
analysis on Fe lines, Israelian et al. (2004) derive a higher grav- 
ity (logg = 1.5), thus indicating that LTE analyses underesti- 
mate the gravity by about 1 dex for this star. 

Stars that may also be located above the tip of the red- 
giant branch (predicted to occur at logg = 0.9 in a 0.8 M , 
[Fe/H]= -3.8 model, as indicated by the horizontal mark on 
the evolutionary tr ack of Fig. [5} also include CS 30314-067 
dAokietalJl207)2cT) . CS 22948-027 dPreston & Snedenl 1200 ll 
but here again Barbuy et al 2005, obtain a l arger log g value 
of 18), and CS 30301-015 (lAoki et al 12002(1 . 

Mau ron et all J2004E005I) list about 120 cool carbon stars 
in the halo (with JHK colors very similar to those of CS 30322- 
023), located at quite large galactocentric distances (up to 
120 kpc!). They were not plotted in Fig.0 however, as neither 
their metallicities nor their gravities are known. 
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Table 4. Distance of CS 30322-023 as a function of the adopted mass, assuming \ogg = -0.3, r ei f = 4100 K and m\, \ = 11.7 
(corresponding to BCy = -0.3 and A v = 0.18 mag). Solutions corresponding to this set of parameters are listed under the header 
'set (1)', whereas an alternative solution, accounting for uncertainties in gravity and temperature, is listed under the header 'set 
(2)' (log g = 0.7, r eff = 4350 K and m bol = 1 1 .7). The ages up to the TP-AGB phase are read from Fig.|H]for the 0.8, 4, and 9 M Q 
models. 



M/M 


log L/L 




d (kpc) 


age (y) 




set (1) 


(2) 


(1) 


(2) 


(1) 


(2) 




10 


5.15 


4.24 


-8.1 


-5.9 


92 


33 




9 


5.10 


4.20 


-8.0 


-5.8 


87 


31 


2.5 10 7 


5 


4.84 


3.94 


-7.4 


-5.1 


65 


23 




4 


4.75 


3.85 


-7.1 


-4.9 


58 


21 


1.6 10 s 


3 


4.62 


3.72 


-6.8 


-4.6 


50 


18 




1 


4.15 


3.24 


-5.6 


-3.4 


29 


10 




0.8 


4.05 


3.14 


-5.4 


-3.1 


26 


9 


13.8 10" 



The low gravity derived for CS 30322-023 appears, at first 
sight, quite puzzling. A comparison with low-metallicity stel- 
lar evolution tracks reveals, however, that both the gravity and 
effective temperature of CS 30322-023 are in fact typical of 
low-metallicity red-giant stars. Fig.0]displays the evolutionary 
tracks in a gravity - effective temperature plane for a 0.8 M Q 
model of metallicity [Fe/H]= -3.8, and for intermediate-mass 
stars of 4, 7, 9 and 12 M and Z = 1 0~ 5 ([Fe/H]= -3.3), com- 
puted with STAREVOL JSiesifeOOd) all the way from the pre- 
main sequence to the thermally-pulsing AGB stage (for the 0.8, 
4, 7 and 9 M models), and to the ignition of n eon for the 12M 
model. For these computations, we use the iReimersl Jl975t) 
mass- loss rate up to the end of cor e He burning, and then switch 
to the Vassili adis & Woodl dl993h prescription during the AGB 
phase. The region around \ogg - and logoff = 3.61 is 
reached by all five models, thus providing a non-unique so- 
lution! However, each of these possible solutions presents its 
own difficulties. In the case of the 0.8 M model, the evolu- 
tionary timescales may be a problem. Although the model star 
has a total lifetime of 13.8 Gyr (about the age of the Galaxy), 
the time spent by such stars in the TP-AGB phase might seem 
prohibitively short to ensure that some are actually observed. 
We will show in Sect.|5]that this is not the case, since the short 
time span on the TP-AGB is more than offset by the large vol- 
ume within which these very luminous stars may be detected. 
In the case of the intermediate or massive stars, one may won- 
der whether such massive short-lived stars could still be found 
in the galactic halo, especially since the distance estimates for 
CS 30322-023 presented in Table@]places it at the outskirts of 
the halo. 

4.2. Observational evidence 

The variability observed for the Ha line profile suggests the 
presence of a stellar wind (Fig.|6j, and provides a further argu- 
ment in favour of the very evolved nature of CS 30322-023. 

Unfortunately, there is no available information on possible 
optical variability, which should be present if CS 30322-023 
lies high on the AGB. The profile resulting from the correla- 
tion of the spectrum with the Arcturus template used to de- 
rive the radial velocity also does not exhibit the asymmetric 
shape expected for Mira variables experiencing shock waves 
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Fig. 4. Evolutionary tracks (in a logg - T e g plane) of a 0.8 M 
model at [Fe/H] = -3.8 (black solid line) and of a 4 M Q (cyan 
short-dashed line), 7 M (red long-dashed line), 9 M (green 
dot-dashed), and 12 M (blue solid line), at met allicity Z = 
10~ 5 ([Fe/H]^ -3.3), computed with STAREVOL llSiessl2006h 
from the main sequence to the thermally-pulsing AGB phase 
(for the 0.8, 4, 7 and 9 M Q models), and to neon ignition for the 
12 M model. The position of CS 30322-023 is represented by 
the asterisk. 

in their photospheres dAlvarez et al.l200lh . Nonetheless, as we 
now show, the available radial-velocity data (Tabled exhibit 
properties that appear to be typical of low-metallicity pulsating 
stars. 

It was shown in Sect. |2] that the radial velocity standard 
deviation (1.24 km s ) is considerably larger than the average 
uncertainty, amounting to 0. 1 8 km s _1 . These variations may be 
caused either by the pulsation of the stellar envelope or by an 
orbital motion in the presence of a companion (or a combina- 
tion of both). Let us first stress that the stellar properties listed 
in Table 0] put stringent constraints on the admissible orbital 
periods. For an 0.8 M Q star, the stellar radius corresponds to 
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Fig. 5. Same as Fig.|4]for the 0.8 M Q model in the usual Mm 
- r e ff plane. 
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Fig. 6. Sample spectra around the Ha line (the flux is in ar- 
bitrary units), obtained on 2002/04/24 and 2002/05/11. The 
change in the profile over 17 days reflects some kind of activity, 
such as, for instance, a strong wind. 



210 Rq with atmospheric parameter set (1) [i.e., T e g = 4100 K 
and logg = -0.3], and to 65 R Q for set (2) [i.e., T eB = 4350 K 
and \ogg = 0.7]. In a binary system involving two stars of 
0.8 M Q (the primary being the observed low-metallicity star 
assumed to be a long-lived low-mass star and the secondary a 
white dwarf), the minimum admissible orbital periods corre- 
sponding to CS 30322-023 not filling its Roche lobe amount to 
1180 d and 200 d, respectively. These periods represent abso- 
lute lower limits, since higher-mass stars have larger radii, and 
hence longer threshold orbital periods. 

Interestingly, a sine-wave of period 192 d very nicely fits 
the radial-velocity data at JD > 2453 000 (upper panel of 
Fig. Q. From the above considerations, this signal has too short 
a period to represent an orbital motion; moreover, it does not 
fit the earlier data points. It is likely, therefore, that it repre- 
sents pulsation of the stellar envelope. The pulsation periods of 



low-metallicity stars are shorter than those of solar-metallicity 
Miras, and thei r range matches the 192 d period observed for 
CS 30322-023 JAlvarez et alJI 19971) . The poor fit to the earlier 
data points may either be caused by a sudden phase shift of 
the (possibly semi-regular) pulsation, or by a long-term trend 
associated with a binary motion. 

We now consider the latter possibility. The dashed line in 
the lower panel of Fig. corresponds to a shift of 1.0 km s _I 
over 1000 d. It provides a better fit to the earlier data points; 
such a trend is not unexpected for binary systems involving 
chemically-peculiar red giants such as CS 30322-023. To fix 
the ideas, let us assume that CS 30322-023 is a binary with a 
period of 5200 d, as is typical for barium stars J.Torissenl20 03). 
In this case, our measurements have only sampled a quarter of 
the orbital cycle. In the most favourable case, the observed ve- 
locity range is then equal to the velocity semi-amplitude, K, 
which would amount to 7.2 km s for a circular orbit seen 
edge-on and involving two stars of 0.8 M G (see above). The 
possible long-term trend of 1.3 km s _1 over 1300 d is not in- 
compatible with the above estimate of K, provided that the or- 
bital plane is inclined on the sky. 

To conclude this discussion, there is no firm support in 
favour of the binary nature of CS 30322-023, although in this 
case the absence of evidence is not evidence of absence, espe- 
cially since any orbital motion will be very difficult to disen- 
tangle from the large intrinsic jitter. The binary hypothesis will 
therefore be left open in the discussion of Sect. 

5. Probability of observing an extremely 
metal-poor TP-AGB star 

Since it is very unlikely that a massive or intermediate-mass 
star has been formed in a recent star-forming episode at such 
large galactocentric distances (> 20 kpc, Tabled in the halo of 
the Galaxy, in this section we evaluate the likelihood of an al- 
ternative scenario, that CS 30322-023 is a TP-AGB star formed 
shortly after the birth of the Galaxy. 

We now estimate the relative fraction of halo TP-AGB stars 
in a magnitude-limited sample (such as the HK survey), as 
compared to main-sequence turnoff or horizontal-branch stars 
from a co-eval sample. The assumption of co-evality is proba- 
bly valid as long as one restricts the sample to extremely metal- 
poor stars with [Fe/H] < -3.0. We further assume that a star of 
O.8M with [Fe/H]= -3.8 has a lifetime (up to the planetary 
nebula phase) precisely equal to the age of the galactic halo 
(Sectgll. 1 

In an ideal co-eval sample {i.e., not subject to observational 
selection biases) consisting of the oldest stars in the Galaxy 
with M > M\ ovl , Mi ow being the low-mass cutoff of the pop- 
ulation (to be fixed below), the fraction, F, of TP-AGB stars 
expected in such a sample is: 

F = . ( 4 ) 

1 We assume in this section that the lifetimes derived from the stel- 
lar models with [Fe/H]= -3.8 apply to CS 30322-023, despite the 
small difference in metallicity. 
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tering these various evolutionary phases. These masses may be 
defined from evolutionary time-scale considerations, since all 
of these stars are supposed to be co-eval with the age of the 
Galaxy TQaiaxy Hence, 

roaiaxy = t MS (0.8) + r TO (0.8) + r RGB (0.8) + r HB (0.8) + t AG b(0.8) 
= t M s(M A gb) + tto(M A gb) + trgb(Magb) + Thb(M a gb) 
= t ms (Mh B ) + t t0 (M ub ) + T RCB (M HB ) 
= t ms (M rcb ) + t t0 (M R gb) 

= T MS (M TO ) 

The last equation of the above set indicates that the stars 
now at the main-sequence TO have a mass Mto, so that their 
lifetime on the main sequence equals the age of the Galaxy. The 
mass, Mrgb, of stars now entering the red-giant-branch stage 
is then given by combining the last two equations in the above 
set: 



tms(Mto) - tms(Mrgb) = t to (M rcb ), 



(5) 
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Fig. 7. Radial velocity as a function of Julian date. Top panel: 
A sine- wave of period 191.7 d, fitted on the data points with 
JD > 2 453 000. Lower panel: The same sine-wave superim- 
posed on a linear trend of 1.0 km s _1 per 1000 d, possibly rep- 
resenting an orbital motion. 



where ^imf is the initial mass function. Here we adopt 
^tmf(M) oc AT 1 - 3 for the mass range [0.1, 0.5] M n , and 
TW(M) oc AT 2 3 for the mass range ]0.5, 0.8] M jjCroupa 
2001). In Eq. |4] Magb corresponds to the mass of a halo 
star currently entering the TP-AGB phase. Before the TP- 
AGB phase, the star will have gone successively through the 
main-sequence, the main sequence turn-off (TO), the red-giant- 
branch phase (RGB), the core helium-burning or horizontal- 
branch phase (HB), and the early- AGB phase, the latter two be- 
ing considered as a whole in what follows. Similarly to Magb, 
one may define the threshold masses Afro, Mrgb, and Mhb 
for stars formed at the birth of the Galaxy and currently en- 



or, in other words, the difference cVtms(Mrgb) of the main- 
sequence lifetimes of the stars of masses Mpo and Mrgb is 
just equal to the time spent on the TO by the star of mass Mrgb 
currently entering the RGB. This difference, 5tms(Mrgb), re- 
lates to the mass difference, 5M = Mro - Mrgb, through the 
lifetime, tms, of main-sequence stars: 



T MS OC M 

leading to 



-2.5 



= -2.5 



5M 



TMS M 

Solving for Mrgb then yields: 



Mrgb = M T0 1+0.4 



T t q(Mrgb) 

tms(M to ) 

The other masses are derived in a similar way: 



(6) 



(7) 



(8) 



M HB = Mjo 



Magb = M T0 (1+0.4 



tms(M T o) 
tto(Magb) 



1 + o \ Tt °( Mhb ) + o \ t kgb( M hb) 
tms(M to ) 
trgb(Magb) 
tms(Mto) 



tms(M T o) 



+ 0.4 



|0/| Thb(M AG b) ) 
tms(M to ) 



0.8 = M T0 (1+0.4 



tto(0.8) 
tms(Mto) 



+ 0.4 



trgb(0.8) 
tms(Mto) 



+o.4 Tm(0 : 8) ,+o.4 TAGB , (0 - 8) ) 



tms(Mto) 



tms(M to ) 



The last equation above is used to derive Mro- An iter- 
ative procedure is required, since the lifetime ratios must be 
normalized by the main-sequence lifetime of a star of yet- 
unknown mass, M T o- In the first iteration, the various lifetime 
ratios normalized by the main-sequence lifetime of the 0.8 M 
star (namely 12.2 10 9 y, up to the turn-off at M(, i ~ 3.5) 
are read from Fig. |5] namely t H b(0.8)/t M s(0.8) = 0.006, 
trgb(0.8)/t MS (0.8) = 0.045, and t to (0.8)/tms(0.8) = 0.08. 
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The lifetime spent on the TP- AGB, Tagb> is ver Y sensitive 
to the adopted mass-loss rate prescription. Very little is known 
about mass loss in extremely low-metallicity stars Jvan Loonl 
120061) . I f it is driven by radiation pressure, then ste llar wind 
models ( Kud ritzki et alii 19871 Il99ll Ivan LoonlEoOoT) indicate 
that the rate depends roughly on the square root of the metal- 
licity. Observational support for the sens itivity of mass-loss 
rate to metallicity has been provided by Gro enewegen et all 
ill 9951) . who compare the dust content of LMC and galac- 
tic AGB stars of s i milar pulsational periods. Alternatively, 
IVassiliadis & Woodl i 19931) derived an empirical formulation 
of the mass-loss rate based on the pulsational properties of 
the Mira variables and pulsating OH/IR stars. Their expres- 
sion is not explicitly dependent on metallicity, but it does in- 
clude the effect of the chemical composition through the stellar 
parameters entering the computation of the pulsational period. 
When the mass-loss rate depends upon metallicity, it remains 
quite low (and hence the AGB lifetime quite long) unless third 
dredge-up (3DUP) episodes increase the metal content of the 
envelope and boost the mass-loss rate. This prescription im- 
plies a longer lifetime on the TP-AGB than the Vassiliadis & 
Wood prescription. Therefore, the TP-AGB lifetime displayed 
in Figs. [H] and |9] computed with the Vassiliadis & Wood pre- 
scription, is likely to represent a lower limit on the actual life- 
time. Moreover, in our calculations, only 26 pulses have taken 
place, and the envelope mass still amounts to 0.1 M with a 
mass loss rate of 4 10" 8 M y -1 . The TP-AGB lifetime dis- 
played in Fig. [8] for the 0.8 M Q model therefore represents a 
lower bound to the actual lifetime, which may be expected to 
be longer by a few 10 6 y. 

Given the uncertainty affecting tagb, it will be kept as a 
free parameter for the remainder of this section. Nevertheless, 
as an illustrative example, conservatively adopting tagb(0.8) = 
10 7 y yields t A gb(0.8)/t M s(0.8) = 8 10~ 4 , hence M T0 = 
0.760 Mq. From Eq.0 such a star has a main-sequence lifetime 
1.5 10 9 y longer than the 0.8 M one, or t M s(A^to)/tms(0.8) = 
1.12. This factor is thus applied to the above lifetime ratios 
to normalize them by tms(A^to), as required. The adopted ra- 
tios are listed in Table [5] This second (and last) iteration then 
yields Mjq = 0.764 M . Since all five threshold masses are 
very close to one another, the normalized lifetimes appearing 
in the equation set above will all be taken equal to the values 
corresponding to the 0.8 M model listed in Table|5] 

To describe a real stellar sample (which in this case is 
magnitude limited), the observational censorship needs to be 
included, implying that the sampling volume will depend on 
the luminosity of the various classes of stars, thus largely 
favouring bright stars, such as AGB stars, over main-sequence 
stars. In the following, the integral appearing in the denom- 
inator of Eq. @] will be split into five different components 
- thermally-pulsing AGB (TP-AGB) stars, horizontal-branch 
and early AGB stars (collectively referred to as HB), main- 
sequence turnoff stars (TO), and finally main-sequence stars 
(MS). The lower limit of the main-sequence mass range con- 
sidered in Eq.0]is denoted M\ ow ; however, as shown below, it 
has no significant impact on the final result. 



it- 

ft v 




Fig. 8. Same as Fig.@]for the ages of the five models as a func- 
tion of the effective gravity. The age is counted backward from 
the last computed model. The occurrence of thermal pulses is 
evident from the rapid oscillations of log g near the end of the 
evolution. The large dots on the 0.8 M track mark the tran- 
sition between different evolutionary phases: main sequence 
(labelled MS), red giant branch (RGB), core helium-burning 
(He-B), early AGB (E-AGB), and thermally-pulsing AGB (TP- 
AGB). 

With these conventions, Eq.|4]can be rewritten 

r 0-&M B 



Pacb 



P 6 r M '+ 

- (= j 72 vol,/ 



,i J M M Vimf{M) AM 



(9) 



where n m \ is the number of halo stars in the volume sampled 
for the various subtypes of stars (namely TP-AGB, HB, RGB, 
TO, and MS; M x = M low ,M 2 = M TO ..., M 6 = 0.8 M Q ), given 
the fact that the HK sample is limited to the 11 < B < 16 mag- 
nitude range (Bee rs et alJl992l) . The adopted volume is a cone 
along the line of sight to CS 30322-023, truncated at a mini- 
mum distance fixed by B = 11, and the absolute B magnitude 
of the considered stellar subtype, and at a maximum distance 
similarly fixed by the condition B — 16. The halo star number 
density is assumed to decrease as [R 2 + (z/0.76) 2 ]~ 244/2 , R be- 
ing the galactocentric distance measured in the galactic plane 
and z the distance from the plane, according to the model of 
iRohinet alJd2003h . 

The distance cut-offs, flf m m,max, f° r the RGB, HB, and AGB 
stars are listed in Table [5] along with the adopted bolometric 
magnitudes and effective temperatures (from the STAREVOL 
tracks discussed in Sect.0]and displayed in Fig.|5}, bolometric 
corrections and total galactic absorption in the B band along 
the line of sight to CS 30322-023, according to the discussion 
in Sect. |3 For TO and MS stars, galactic absorption is assumed 
to be negligible, since the volume sampled is small. The bolo- 
metric correction in the B band is taken from a grid of MARCS 
models with [Fe/H] = -4, and a range of gravities and effective 
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Fig. 9. Same as Fig.0for the 0.8 M model, but expressed as 
a function of M\, \. The lifetimes in the different evolutionary 
stages are used in Sect. [5] Note that the AGB lifetime provided 
by this figure is a lower limit to the actual value, since only 
26 pulses have been computed, and the remaining envelope 
mass amounts to 0.1 M . A few million years should thus be 
added to the AGB lifetime displayed on the figure. Moreover, 
the Vassiliadis & Wood (1993) prescription for the mass loss 
used in these computations (see text) also shortens the AGB 
lifetime. 



temperatures. The y do not diffe r much from those provided in 
Table 6 of lAlonso etail dl999l) for giants with [Fe/H] = -3. 
The bolometric correction for AGB stars has been taken iden- 
tical to that derived for CS 30322-023 (Sect.0. 

Main-sequence stars are sampled in an exceedingly small 
volume as compared to the other, more luminous classes, as 
may be judged from Table |5] Th ose stars s hould thus be very 
rare among the HK survey jBeers et all 19851) since the volume 
effect largely offsets the effect of the increasing IMF, hence 
the first term (i = 1) of the sum in Eq. |9]is neglected in the 
following. 

The number of stars in a truncated cone along the line of 
sight to CS 30322-023, with an opening angle of a (this pa- 
rameter, appearing in both the numerator and denominator of 
Eq.|9] has no impact on the final result) thus becomes 

« vo i(M) =7rtg 2 apo [R 2 + (z/0.76) 2 ) ' d 2 dd (10) 

where R = ro - c/cos b (or d cos b - ro, whichever is positive) is 
the galactocentric distance measured in the galactic plane, z = 
d sinb the distance from the plane, b = -46.3° is the galactic 
latitude of CS 30322-023 (its galactic longitude, / = -6.8°, is 
sufficiently close to to be neglected in the above formula), 
ro ~ 8 kpc is the galactocentric distance of the Sun, and po is 
the central galactic density of extremely metal-poor halo stars. 




Fig. 10. Fraction of AGB stars (solid line), HB stars (dotted 
line), RGB stars (short-dashed line), and TO stars (long-dashed 
line) in a co-eval sample born in the early Galaxy. 



The latter parameter will later cancel, and thus has no impact 
on the final result. 

The integrals of Eq. [TO] have been computed numeri- 
cally; the results are listed in Table [5] The final proportion 
of AGB stars in the co-eval sample is displayed in Fig. ^3 
Given the fact that a reasonable estimate for the fractional TP- 
AGB lifetime, tagb(0.8)/t M s(M T o), is on the order of 10~ 3 (= 
10 7 /10 10 y/y), the main conclusion from the present analysis 
is that it is not at all unlikely to observe an extr emely low- 
metall icity TP-AGB star in a co-eval sample. The B eers et all 
( 1992) sample is most likely not a co-eval sample, as it includes 
stars with metallicities larger than that of CS 30322-023 , which 
must have formed somewhat more recently. In a sample con- 
taining stars formed over a long period of time (with a constant 
star formation rate), the relative fractions of stars of various 
kinds are proportional to their lifetimes. Since the lifetime of 
stars at the main-sequence turnoff is about 7 times longer than 
the AGB lifetime, it is not surprising that Bee rs et al T ill 9921) 
sample contains a much larger fraction of TO stars than pre- 
dicted under the co-evality hypothesis. Nevertheless, it should 
be stres sed that the above analysis was not designed to mim- 
ick the iBeers et alJ £l992) sample (this would require a de- 
tailed knowledge of the star formation history in the halo, 
which is beyond the scope of the present paper), but rather to 
roughly evaluate the likelihood of finding an extremely low- 
metallicity TP-AGB star. Incidentally, an analysis similar to 
the on e presented here, based o n the Bahcall-Son eira Galaxy 
mode l ( Bahcal l & Soneirall980l) . was presented bv lBeers et alJ 
(1985), who arrived at the same conclusion - a magnitude- 
limited sample of halo stars to apparent magnitude B — 16 is, 
by far, dominated by giant stars. 
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Table 5. Sampling distances of the HK survey for the various stellar subtypes and fractional lifetimes t(0.8)/t M s(M T o) for 
the various evolutionary phases of the 0.8 M G model of [Fe/H]= -3.8. The lifetime in the AGB phase is subject to various 
uncertainties, and has therefore not been listed. BCg and A B are, respectively, the bolometric correction and total absorption in 
the B band. 



type 


t,(0.8)/t ms (M to ) 


log T cff 


M bo i 


BC B 


A B 


^min 


^max 


w vol/ w vol,TO 


i 












(kpc) 


(kpc) 




AGB 




3.61 


-5 


-1.6 


0.2 


6.91 


69.1 


274 


HB + E-AGB 


5.6 10- 3 


3.85 


0.5 


-0.5 


0.2 


0.91 


9.12 


40.8 


RGB 


4.0 10~ 2 


3.71 


0.5 


-1.0 


0.2 


0.72 


7.24 


24.7 


TO 


6.9 10~ 2 


3.85 


3.5 


-0.5 


0.0 


0.25 


2.51 


1 


MS 


0.88 


3.81 


5 


-0.5 


0.0 


0.12 


1.23 





Table 6. Abundances for CS 30322-023 obtained with the 
plane-parallel model. N is the number of lines used in the anal- 
ysis; 'syn' indicates that spectral synthesis has been used to 
derive the abundance. The column log e(X) lists the abundance 
in the scale where log e(H) = 12. Ab undances are normalized 
by the solar data from Asclun det alJ (|2P05). 



Element 


N 


loge(X) 




[X/H] 


[X/Fe] 


C(CH) 


(syn) 


5.60 


_ 


-2.79 


0.60 


N(CN) 


(syn) 


7.20 


_ 


-0.58 


2.81 


OI 


l(syn) 


5.90 


_ 


-2.76 


0.63 


Nal 


4 


4.07 


0.24 


-2.10 


1.29 


Mgl 


4 


4.94 


0.10 


-2.59 


0.80 


All 


2 


3.70 


0.3 


-2.67 


0.72 


Si I 


1 


4.70 


0.2 


-2.81 


0.58 


KI 


l(syn) 


<3.00 




-2.08 


1.31 


Cal 


10 


3.37 


0.12 


-2.94 


0.45 


Sell 


4 


-0.12 


0.17 


-3.17 


0.24 


Til 


10 


1.31 


0.19 


-3.59 


-0.20 


Till 


14 


1.72 


0.13 


-3.18 


0.23 


VII 


3 (syn) 


0.40 


0.1 


-3.60 


-0.19 


CrI 


9 


1.66 


0.22 


-3.98 


-0.59 


Mnl 


6 


1.66 


0.28 


-3.73 


-0.34 


Fel 


66 


4.06 


0.18 


-3.39 





Fell 


7 


4.04 


0.19 


-3.41 





Col 


1 


1.33 




-3.59 


-0.20 


Nil 


6 


2.56 


0.18 


-3.67 


-0.28 


Cul 


l(syn) 


-0.10 




-4.31 


-0.92 


Znl 


2 


1.23 


0.04 


-3.37 


0.02 


SrII 


l(syn) 


-1.00 




-3.92 


-0.51 


YII 


6 


-1.56 


0.15 


-3.77 


-0.36 


Zrll 


5 


-0.68 


0.30 


-3.27 


0.14 


Tel 


l(syn) 


<-1.40 








Rul 


l(syn) 


<-0.50 




-2.34 


1.05 


RhI 


l(syn) 


<-2.00 




-3.12 


0.27 


Ball 


3 (syn) 


-0.70 


0.1 


-2.87 


0.54 


La II 


8(syn) 


-1.80 


0.1 


-2.93 


0.48 


Ce II 


18 


-1.22 


0.24 


-2.80 


0.61 


Pr II 


3(syn) 


-2.50 


0.10 


-3.21 


0.20 


Ndll 


24 


-1.34 


0.29 


-2.79 


0.62 


Smll 


4(syn) 


-2.30 


0.2 


-3.31 


0.10 


Eu II 


l(syn) 


-3.40 




-3.92 


-0.51 


Gdll 


l(syn) 


<-2.30 




-3.42 


-0.01 


Tbll 


l(syn) 


<-3.20 




-3.48 


-0.07 


Dy II 


2(syn) 


-2.40 


0.1 


-3.54 


-0.13 


ErII 


l(syn) 


-2.56 




-3.49 


-0.08 


Tmll 


l(syn) 


<-3.50 




-3.50 


-0.09 


Ybll 


1 


-2.10 




-3.18 


0.23 


Hfll 


2(syn) 


-2.00 


0.1 


-2.88 


0.53 


WI 


l(syn) 


<-1.50 




-2.61 


0.78 


Pbl 


2(syn) 


0.10 


0.2 


-1.90 


1.49 



l2 C/ 13 C = 4 ± 1 




[FVH] 



Fig. 11. [C/H] as a function of metallicity for metal-poor stars, 
collected from the literature (see Masseron 2006, for detailed 
references). The lower diagonal sequence corresponds to the 
trend expected from galactic chemical evolu tion for unpollute d 
stars. Two of the stars from the sample ol'iMasseron! |2006), 
among which are CS 30322-023 (asterisk), fall on this galac- 
tic sequence. The picture changes, however, when considering 
C+N (Fig.IH. 



6. Abundances 

Because of the general crowding of the spectra for cool, 
evolved stars with molecular lines, the analysis of stars as 
CS 30322-023 is very challenging. Thanks to the low metal- 
licity, however, abundances could be obtained for 35 species, 
while seven received upper limits on their abundances with 
the plane-parallel (PP) model. These results are summarised 
in Table [6] Hyper-fine structure has been included for Ba, La, 
Pr, Eu, Tb, Dy, Er, Lu, and Pb. C and N abundances are derived 
from the CH and CN lines. The O abundance is derived from 
the 630.0 nm [O I] forbidden line which, among the O abun- 
dance diagnostics, is expected to be the least affected one by 
non-LTE conditions ( Israelian et al. 2004), although it is sensi- 
tive to 3D granulation effects ( Shchuki na et al.l2 005). 

Although the gravity has been derived from the Fe I/Fe II 
ionisation balance, a large discrepancy (about 0.4 dex) is nev- 
ertheless found between the abundances derived from the Ti I 
and Ti II lines. This discrepancy might be ascribed to vari- 
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Fig. 12. Same as Fig.^2f° r [N/H] as a function of metallicity 
for CEMP stars collected from the literature (see Sect.|4]for the 
full reference list). As usual, CS 30322-023 is represented by 
an asterisk. 
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Fig. 14. [(C+N)/H] as a function of metallicity for the same 
sample of CEMP stars as in Fig. 1121 
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Fig. 13. [O/Fe] as a function of [Fe/HJ for CS 30322- 
023 (a sterisk), CS 29498- 043 lAoki et aljEool) . HE 0107- 
5240 (Besselletal. 2004) and for unevolved dwarf stars 
from Israelian et al. (2001) ( open squar es). ICavrel et al. (2004) 
(filled circles'). and lNissen et alJ J2002I) (filled diamonds). The 
filled symbols correspond to [O/Fe] values derived from the [O 
I] line, while open symbols correspond to oxygen abundances 
derived from the triplet or OH li nes. No 3D correc tions have 
been applied. Figure adapted from lAoki et al.ld2004l) . 




Fig. 15. Abundance pattern of CS 30322-23 compared to two 
model predictions for the s- and r-processes. The s-process 
abundance pattern has been computed by post-processing an 
AGB model o f mass 0.8 M H and metallic ity [Fe/H] = -3.8, 
as described in iGorielv & Mowlavil feOOCl) . The r-process pat- 
tern corresponds to th e scaled solar r-process abundances, from 
iPaime&Beerl ill 9941) . 



ous causes: (i) the use of a plane-parallel instead of a spher- 
ical mo del (but see the discus s ion in Sect. Et, (ii) NLTE 
effects JThevenin & Idiarl ll999Hlsraelian et alTl200ll 12004 
iKorn et alJ2003h . (iiil 3D~effects llAsplundl2005l) . or (iv) a sys- 
tematic shift between the scales on which the log gf values for 
Ti I and Ti II are derived. A similar proble m has bee n reported 
for Cr in other cool metal-poor stars Jjohnsonl2002l) . The first 
three possibilities should be investigated in the future by using 
tailored 3D, NLTE models. 

The uncertainty quoted in Table [6] in the column labelled 
cr does not include such possible systematic errors, as it only 
refers to the standard deviation of the abundances derived from 
the different lines used for a given species. When only one line 
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is available, spectrum synthesis has been used, and a 0. 1 dex 
uncertainty should in this case be considered as typical, arising 
from errors associated with the location of the continuum and 
on the line-fitting process. 

Table [6]reveals that, quite remarkably, the C enhancement 
in CS 30322-023 is only moderate, in contrast to the N en- 
hancement (Figs.fTT1and ll2> . Na is strongly enhanced as well. 
The C/O ratio is only 0.50. The C and O abundances are in 
fact compatible with the trend due to the chemical evolution 
of the Galaxy, as displayed in Fig.^2f° r C and in Fig.ll3lfor 
O. However, because the carbon isotopic ratio is low ( 12 C/ 13 C 
= 4), it is unlikely that C has retained its initial value in the 
envelope of CS 30322-023. As further discussed in Sect. 
the combination of a low 12 C/ 13 C ratio and large N overabun- 
dance strongly suggests that pollution of the atmosphere by 
material processed by the CN cycle has occurred. In this re- 
spect, it is worth noting that CS 30322-023 falls among the 
other CEMP stars in a diagram showing [(C+N)/H] as a func- 
tion of metallicity (Fig.[H). The a-elements in CS 30322-023 
do not exhibit any overabundances with respect to the values 
expected for its metallicity from galactic chemical evolution. 
There is also no observed enhancement of r-process elements, 
whereas the s-process elements belonging to the second and 
third peaks are clearly enhanced (Fig. I15i. Tc I and Tc II lines 
have been searched for, and not found. As shown on Fig.^] an 
upper limit log e(Tc) < -1.4 (translating to e(Tc)/e(Zr) < 0.2) 
may be derived from the Tc I /1426.23 nm line. This upper 
limit on the Tc abundance is not meaningful, however, since 
e(Tc)/e(Zr) ~ 0.001 is expected from s-proces s nucleosynthe- 
sis (see Fig. 15 of iGorielv & MowlavilEoOoT) . The difficulty 
of detecting the Tc I line in this star results from its rather 
high temperature (as compared to solar-metallicity AGB stars, 
which are much cooler), implying that Tc is mostly ionized. 
Unfortunately, the strongest Tc II lines lie in the far ultravio- 
let; perhaps these lines could be searched for from space in the 
future. 

The abundance pattern for heavy elements is presented 
in Fig. [21 along with the scaled solar r-process abundance 
pattern, and an s-process abundance pattern expected from 
proton-mixing, forming a 13 C pocket below the He-H discon- 
tinuity in a 0.8 M n AGB star of metallicity [Fe/H] = -3.8 
( IGorielv & MowlavilEoOOl) . With the exception of Sr and Y, 
all observed abundances almost perfectly match the predictions 
of a pure s-process operating at low metallicit y in the frame- 
work of the proton-mix ing scenario dStraniero et alJ 1 19951: 
IGorielv & Mowlavil 2000). The s-process abundance pattern in 
CS 30322-023 is at the lowest metallicity observed to date. 
As expected, it differs from its so lar-metallicity co unterpart 
due to the large Pb overabundance JGorielv & Mowlavil f2000 ; 
iBusso et alJboOlt IGorielv" & Siess 2001), which matches the 
value observed in CS 30322-023 (Figs.fT31andlT7t. th us making 
this star a ne w member of the class of lead stars ( Va n Eck et alJ 
|2001| ]2003). The large spread in the Pb abundances observed 
among CEMP stars (Fig. 1171 is surprising, however, as is fur- 
ther discussed in Sect. [7] The parameter responsible for this 
scatter remains unidentified, as the scatter does not correlate 
with either [C/Fe] or N/(C+N). 
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Fig. 16. The observed spectrum of CS 30322-023 (crosses) 
around the Tc I /I426.23 nm line, with three synthetic spectra: 
without Tc (blue solid line), with log e(Tc) = -1 .4 (red dashed 
line), and with log e(Tc) = -1.0 (cyan dot-dashed line). 
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Fig. 17. The [Pb/Ba] ratio as a function of metallicity for the 
sample of CEMP stars (excluding r-process-enhanced stars) 
taken from the literature (see the references listed in Sect.|4j. 
As before, CS 30322-023 is represented by an asterisk. 



To estimate the global overabundance of s-process elements 
in CEMP stars - when evaluating, for instance, its possible cor- 
relation with the C and N abundances - it is important to con- 
sider not only overabundances of elements (such as Ba) belong- 
ing to the second s-process peak, but also overabundances of el- 
ements (such as Pb) from the third peak. Failure to do so would 
largely underestimate the s-process overabundance whenever 
lead stars are involved. This is very clear from Fig.^] Fig. 1181 
reveals that there is indeed a general tendency for ([Ba/Fe] + 
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Fig. 18. The global overabundance of s-process elements in 
CEMP stars (excluding r-process-enhanced stars), expressed 
in terms of [Ba/Fe] + [Pb/Fe], as a function of [(C+N)/H]. 
Although the general trend may be explained by the dilution 
of the processed material in the envelope, there remains a sig- 
nificant - unexplained - scatter in the s-process abundances (at 
a given [(C+N)/H]). 
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Fig. 19. Same as Fig. ^1 f° r [Pb/Ba] as a function of 
[(C+N)/H]. 

[Pb/Fe]) to increase with [(C+N)/H] (this trend may in fact be 
ascribed to the dilution of the processed material in a more 
or less massive envelope), although there remains a consider- 
able scatter. On the contrary, [(C+N)/H] exhibits little correla- 
tion with the s-process efficiency expressed in terms of [Pb/Ba] 

(Fig.ua. 

In absolute terms, the s-process overabundance ([Ba/Fe] + 
[Pb/Fe]) in CS 30322-023 is thus moderate, when compared 



to other CEMP stars (Fig. 1 1 81 . which is suggestive of a large 
dilution factor of the processed material in the envelope. 

Considering the sum Ba+Pb, instead of solely Ba, is also of 
interest when discussing the [Ba/Eu] ratio. Fig.l20lshows that, 
in a [Ba/Eu] vs. [C/Fe] diagram, there are several stars with 
[Ba/Eu] ratios that are intermediate between the solar s-process 
and solar r- process ratios. These stars have been coined 's+r' 
or 's/r' stars faill et alfcOOOHCohen et all2003aHSivarani et alJ 
l2004lBarbuv et alJ2005Ulvans et alJ2005tlJonsell et all200d) : 
their [Ba/Eu] ratios do indeed seem peculiar when compared 
to the solar s-process ratio (Fig. l20t . Detailed s-process com- 
putations in low-metallicity AGB stars dGorielv & SiesJ200d) 
have shown that the [Ba/Eu] ratio predicted under such circum- 
stances are quite different from the solar s-process Ba/Eu ratio 
(upper panel of Fig.l20>. and matches at least some of the ob- 
served [Ba/Eu] values. When considering the [(Ba+Pb)/Eu] ra- 
tio, instead of the [Ba/Eu] ratio, s-process predictions become 
more robust (in the sense that the different environments yield 
similar values, within 0.6 dex). Such predictions may account 
for some of the stars displayed in the lower panel of Fig. 
but the need for a separate class of 's+r' ('s/r') stars is even 
stronger, as they separate even more clearly from pure s- and 
pure r-stars when considering [(Ba+Pb )/Eu] instead of [Ba/Eu]. 

7. Observed abundances and possible models 

In solar-metallicity AGB stars, C is brought to the surface by 
the third dredge-up, hence the C/O ratio for a given star is 
expected to increase steadily with time. Nitrogen is not en- 
hanced, unless the operation of hot-bottom burning (HBB) in 
intermediate-mass AGB stars has converted C into N. The sit- 
uation in low-metallicity AGB stars is not as simple, as several 
specific mixing processes with potentially large impact on the 
surface abundances have been found to operate in these stars. 
Since CS 30322-023 is believed to be a low-mass, extremely 
low-metallicity TP-AGB star, as we have argued in Sect. |5| 
we now briefly describe the various mixing processes in low- 
metallicity stars that have been specifically discussed in the lit- 
erature: 

- Some deep mixing is believed to occur in low-mass, low- 
metallicity stars on the first giant branch, after the luminos- 
ity bump. This deep mixing brings to the surface hydrogen- 
burning ashes, including Na and Al in s ome instances, and 
is a signature of hot hy drogen burning JCharbonnelll 1 994 
ll995llSnite et all2006h 

- He-flash-driven deep mixing (He-FDDM) occurs at the tip 
of the RGB when the helium-flash co nvective zone ex tends 
into the tail of the H - burning shell JFuiimo to et all 1 1990: 
Hollowell et aljfl99(it iFuiimoto et alJl2000[ ISchlattl et alJ 
200l[l2002Ulwamoto et alJ2004ISuda et alJ2004 . The in- 
jection of protons in the region of He-burning produces a 
H-flash and causes the convective region to split. The lower 
convective shell is powered by He burning while the up- 
per one is driven by H burning. The structural readjustment 
subsequent to the H-flash leads to a CN-strong star, through 
the merging of the H-driven convective shell (formerly pol- 
luted with 3ff- 12 C) with the outer envelope. It is important 
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Fig. 20. Upper panel: [Ba/Eu] as a function of [C/Fe]. Lower 
panel: [(Ba+Pb)/Eu] as a function of [C/Fe]. Triangles: r- 
process-enhanced stars; circles: s-process-enhanced stars; as- 
terisk: CS 30322-023. Open symbols correspond to stars with 
either no Pb data or with only an upper limit available. The 
solar [Ba/Eu] a nd rPb/Eu] rat ios ascribed to the r- and s- 
processes (from Gori elviri999l) are displayed as dotted lines, 
along with the maximum r-proces s and min imum s -process so- 
lar ratios, taken from Table 3 of iGorieivl (1999). Predictions 
from the 0.8 M Q , [Fe/H] = -3.8 model (see Fig.fBl and from 
the s-process o perating in massive AGB stars (from Fig. 6 of 
iGorielv & Sies s 2006) are also shown. 



to note, however, that the atmosphere of such a CN-strong 
star is not expected to be enriched in s-process elements. If 
the injection of protons in the He flash leads to the produc- 
tion of 13 C, and subsequently of s-elements, these heavy 
elements will remain trapped in the He-burning convective 
region, which the envelope never penetrates; 



Very similar to the previous scenario is the mixing of pro- 
tons into the thermal pulse at the beginning of the TP-AGB 
phase, which also leads to a CN-strong star iFuiimoto et alJ 
200Q| ISuda et alJl2004l) . In contrast with the previous sit- 
uation, however, this process might produce s-process el- 
ements, but this needs to be confirmed (Aikawa, priv. 
comm.); 

"Carbon injection" is a mi xing process encountered ex- 
clusiv ely in Z = models dChieffi et alJl200lt ISiess et alJ 
2002). It occurs at the end of the early AGB phase 
when the He-burning shell (He-BS), while approaching the 
hydrogen-burning shell, is subject to small- amplitude in- 
stabilities. In this process, after the development of a weak 
thermal pulse in the He-BS, a secondary convective zone 
forms at the He-H discontinuity. It grows in mass, pene- 
trates downwards into the C-rich layers left over from the 
pulse, and initiates a H-flash by activation of CNO reac- 
tions. The subsequent structural readjustment causes the 
envelope to penetrate in this secondary convective zone, 
leading to a surface CN-enrichment; 
The "classical" proton-mixing scenario assumes that some 
extra mixing takes place at the base of the convective en- 
velope during the third dredge-up. It is responsible for 
the injection of protons into the C-rich layers left by the 
pulse, leading to the formation of a 13 C pocket. The radia- 
tive burning of 13 C during the interpulse phase then leads 
to the production of s-process elements which, after pass- 
ing through the pulse, are brought to the surface by the 
third dredge-up JStraniero et alJll99^ IGorielv & Mowlavl 
2000). 



In addition to these extra-mixing mechanisms, HBB can 
also significantly modify the observed surface abundances, at 
least in intermediate-mass stars, when H-burning occurs at the 
base of the convective envelope, and partly converts 12 C into 
14 N. HBB can also activate the NeNa and MgAl chains, and 
therefore could be responsible for 23 Na and 27 A1 enrichment. 
In Pop III stars, HBB can be activa ted in stars with masses as 
low as M ^ 2M Q dSiess et all2002l) . 

Signatures of these mixing and HBB processes generally 
are a combined enhancement of hydrogen-burning products 
(low 12 C/ 13 C, high N, and sometimes high Na) and helium- 
burning products ( 12 C and possibly s-process elements) in the 
stellar envelope. Since the stellar envelope is enriched by pri- 
mary 12 C, C-enhancement factors of several dex are expected. 
For example, in the 0.8 M model described above, mixing be- 
tween the He-flash and the stellar envelope increases the en- 
velope carbon content to [C/H]^ +1.0. Such large overabun- 
dances are not confirmed by observations thus far (see Fig. ll II . 
The abundances in CS 30322-023 are instead suggestive of CN- 
cycling of formerly C-rich matter, perhaps through HBB. Since 
HBB is not expected to occur in a 0.8 M star, we are thus 
left with the following alternatives. Either some extra mixing 
occurs in low-mass stars that results in H-burning similar to 
HBB, or an enrichment by an intermediate-mass AGB com- 
panion which has undergone HBB has polluted the presently- 
observed star. The latter suggestion must await firm evidence 
for orbital motion, which is lacking at present. 
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Further constraints on the the nucleosynthetic history of 
this star come from the s-process enrichment observed in 
CS 30322-023. The s-process in AGB stars is so far under- 
stood by the partial mixing of protons in C-rich layers, fol- 
lowed by the 3DUP of this material (Strani eroet all 1 19951 
Idorielv & Mowlavil2000h . As shown in Fig.Q3] the CS 30322- 
023 s-process abundance pattern is in excellent agreement with 
the so-called partial-mixing model. However, in such a sce- 
nario, the s-process surface enrichment is inevitably accom- 
panied by a large C enrichment due to the 3DUP of primary 
C made in the thermal pulse, in contradiction to the low C- 
to-s-process overabundance exhibited by CS 30322-023. None 
of the present s-process models could possibly explain such a 
ratio; for this reason CS 30322-023 is quite simi lar to the C- 
poor s -rich post-AGB star V453 Oph observed bv lDeroo et alJ 

More generally, the observed relation between the [Pb/Ba] 
ratio and metallicity shown in Fig. El is not reproduced by 
any current models. The traditional proton-mixing scenario 
jGorielv & Mowlavil Eoool) cannot explain the low [Pb/Ba] 
stars observed at low metallicity, whereas the mixing of protons 
in the thermal pulse at the beginning of t he TP- AGB phas e, ac- 
cordin g to the models of lFuiimoto et alJ ( 120001) and lSuda et alJ 
(2004), predicts a more complicated trend between [Fe/Hl and 
[Pb/Ba], not seen in Fig.^l( com P are with Fig- 7 of ISuda et alJ 

Our new observation , as we ll as the C-poor, s-rich post- 
AGB star of iDeroo et alJ II2005I) . clearly show that our under- 
standing of the mixing mechanism(s) that are responsible, not 
only for the surface enrichment, but also for more specific nu- 
cleosynthesis events such as the s-process, are still far from be- 
ing under control. Our ability to quantitatively explain the ex- 
isting observations remains highly dependent on stellar model- 
ing, and most of all on the description of convective boundaries 
where some extra mixing takes place. A better description of 
the internal transport processes must also be accompanied with 
better knowledge of the mass loss, which indirectly influences 
the thermodynamical conditions at the base of the convective 
envelope and consequently the 3DUP events. 

8. Summary 

The star CS 30322-023 appears to be the most metal-poor 
([Fe/H] = -3.5) s-process-rich star known so far. It is also 
a lead star, with [Pb/Ba] = +0.95. An LTE analysis of the 
Fe I/Fe II ionization balance yields a very low gravity of log g = 
-0.3. Although it is likely that such an LTE analysis underes- 
timates the gravity, a comparison with other CEMP stars ana- 
lyzed with the LTE assumption reveals that CS 30322-023 is 
the most evolved example among known extremely metal-poor 
stars. It is most likely to be a low-mass (~ 0.8 M@) TP-AGB 
star. Larger masses would put CS 30322-023 too far (> 20 kpc) 
in the outskirsts of the halo, where no recent star formation 
is expected. Despite the rather short lifetime of the TP-AGB 
phase, the high luminosity of TP-AGB stars makes it possi- 
ble to sample them in a very large volume. The probability to 
uncover such low-mass, low-metallicity TP-AGB stars is there- 
fore not small, as confirmed by a detailed quantitative analysis. 



We thus expect that additional examples will be revealed in the 
near future, as spectroscopic surveys of the Galaxy push to ever 
fainter limiting magnitudes. 

The observed abundance pattern of CS 30322-023 is, how- 
ever, quite puzzling, in that it exhibits s-process elements with- 
out a strong C overabundance, as would be expected if the 
s-process operates in a He-burning environment. The light- 
element abundance pattern is instead reminiscent of the oper- 
ation of the CN cycle (strong N overabundance, low 12 C/ 13 C 
ratio, Na overabundance). Hot-bottom burning is not expected 
to operate in such low-mass stars. We therefore suggest that 
the enrichment of the envelope with nuclear-burning products 
does not result from the traditional third dredge-up events, but 
rather may be caused by some mixing mechanism that has yet 
to be identified. Although there is no firm evidence at present 
for orbital motion of this star, an alternative possibility is that 
CS 30322-023 resides in a binary system, and has been pol- 
luted by matter from an intermediate-mass companion that pro- 
duced s-process elements and converted C into N during HBB 
episodes. Clearly, full understanding of the abundance pattern 
and apparent evolutionary state of CS 30322-023 will depend 
on future high-resolution studies of stars with similar charac- 
teristics. Abundance determinations duly accounting for NLTE 
and 3D effects, which are suspected to play an important role 
in metal-poor giant stars, should as well bring progress in this 
field. 
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